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We derive effective dynamical theories for metals and BCS superconductors, based on the effective 
action formalism. Both the metallic regime T > Tc and the superconducting regime T -C To are 
studied in the clean and dirty limit. Furthermore, we consider the effect of particle-hole asymmetry 
in the band structure. Using gauge invariance, the electrodynamics of the problem is formulated in 
a transparent way. The effective actions are useful starting points for treating dynamical problems 
involving BCS superconductors. 
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PACS numbers: 74.20.Fg, 71.10.-w 
I. INTRODUCTION 

Dynamical problems in BCS theory! are diverse. They 
include the electromagnetic response of superconduc- 
torsel, relaxation.jphenomena, and collective modes in 
superconduciprsOu, e.g., the Carlson-Goldmanfl and the 
Mooij-SchonO mode. Further examples are the motion 
of topological defects, e.g., vortex motion in bulk sam- 
ples0~£3, quantum tunneling of vorticesta, and thermally 
activated, or quantum phase slips in one dimensional 
wireaL3'E3, as well as fluctuation effects, e.g., corrections 
to the conductivity above TcoE£ni3, the renormalization 
of the critical temperature and of the energy gap of low 
dimensional dirty supcrcond-uctorsEj, and the quantum 
melting of the vortex latticaia. For all these phenomena, 
an effective (simple) theory of weak coupling BCS super- 
conductivity is desirable. However, such an effective the- 
ory is well established only close to and above the critical 
temperature, where Time Dependent Ginzburg-Landau 
(TDGL) theory - although only under severe additional 
restrictions - can be derived from a microscopic starting 
pointEBL Since the. works of Abrahams and ' 
PopovEil, KleinertO, AmbegaokarE 3 ], and Stoo: 
known about the extensions to lower temneraiurc 
controversies concerning the subject persisiO~tJ. 

Previous papers were mostly restricted to the clean 
limit at zero-temperature and neglected the Coulomb in- 
teractionE3~El In the present paper, we extend the ex- 
isting literature in four ways. 1) We fully account for 
the Coulomb interaction between electrons and between 
electrons and the ionic Jellium background. 2) We also 
consider the dirty limit, in which electrons move diffu- 
sively rather than ballistically. 3) We allow for particle- 
hole asymmetry, i.e., the dependence of the density of 
states on energy, as quantified by its derivative at the 
Fermi surface Nq = d e N(ep)- The relevance of particle- 
hole asymmetry for vortex motion and the corresp 
ing Hall-effect has recently been pointed out in Refs. 
4) Finally, a guiding principle is called for, as the ex- 
pansion of the effective action that we will perform (see 
below) is quite involved. We will make extensive use of 
gauge invariance and the corresponding Ward-identities. 
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Although perturbation expansions up to a definite order 
may break gauge-invariance, the Ward identities allow us 
to construct a manifestly gauge-invariant effective action 
within pprt.nrhatinn theory. 

In RefsHjo Galilei invariance in the clean limit was 
used as a guiding principle. We argue that in real su- 
perconductors Galilei invariance can be broken by at 
least four mechanisms. 1) Scattering on impurities and 
phonons break Galilei invariance in a trivial way. 2) 
Galilei invariance would require a perfect quadratic dis- 
persion, which is not realized in the bandstructure of 
usual crystal backgrounds. 3) The coupling to electro- 
magnetism, which is Lorentz invariant, also (weakly) 
breaks Galilei invariance (a vp/c effect). 4) At nonzero 
temperature the gas of normal quasi-particle excitations 
provides a preferred frame of reference which breaks 
Galilei invariance. 

Instead of Galilei invariance, we stress the role of 
gauge invariance. The corresponding Ward identities ex- 
press particle number and charge conservation and can 
be used to rewrite the effective action in a manifestly 
gauge-invariant way. This means that in the supercon- 
ducting phase the effective action depends only on the 
amplitude of the energy gap |A|, the superfluid veloc- 
ity v s = 2^r(V<,5 — — A) and the chemical potential for 
Cooper pairs $ = V — ip/2e. r-Ihis is in accord with 
the Anderson-Higgs mechanismEj: the superconducting 
phase does not appear explicitly in the effective action. 
We will see below that, for example, if the equilibrium 
value of the energy gap |Ao| does not depend on space 
and time coordinates and all fields vary slowly in space 
and time, the action in the particle-hole symmetric case 
has the form 
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Here |A]J is the small deviation of the order parameter 
from the equilibrium value, the x' s are generalized sus- 
ceptibilities, e is the dielectric function of the metal, and 
\x is the magnetic permeability. The terms involving E 
and B include the electronic polarization contributions 
to the electromagnetic fields that will be discussed below 
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in Section III. The term involving the amplitude fluctu- 
ations | Ai| is discussed in Section IV. Most interesting 
are the last two terms: the first describes how the gra- 
dient of the superconducting phase tries to adjust to 
the local vector potential A. The prefactor \l is propor- 
tional to the superfluid density n s and this term is just 
the kinetic energy of the superfluid, related to the DC 
Josephson effect and to London theory. Similarly, the 
less familiar second term describes how the time deriva- 
tive of the superconducting phase tp tries to adjust to the 
local scalar potential V. This is the term that produces 
the AC Josephson effect. The prefactor \J is propor- 
tional to the superfluid density as well. If the critical 
temperature Tc is approached from below, the superfluid 
density vanishes, and the chemical potential for Cooper 
pairs $ and the chemical potential for quasi-particles V 
decouple. This last term also describes how the motion of 
vortices in the mixed state leads to a voltage drop across 
the sample via the time dependence of the phase. 

In case particle-hole symmetry is broken, an additional 
contribution to the action arises, which is proportional 
to the derivative of the density of states at the Fermi 
energy Nq = d e N(eF)- The dimensionless parameter 
that characterizes the amount of particle-hole symmetry 
breaking is 7 = AN( ) /(2\N§), with BCS coupling con- 
stant A. In usual weak coupling superconductors with 
XNq ~ 10 _1 and A/e^ ~ 10~ 3 the parameter 7 is rather 
small, 7 ~ 1CP 2 . Nevertheless, it is important for vor- 
tex motion, and its possible relation to the sign=anomaly 
in the Hall effect has been discussed in Refs.ElEj. The 
general form of the particle-hole asymmetric part in the 
action is 
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(2) 



with r = 7/A0. The physical origin of this term is the 
coupling of the electronic density to the energy gap when 
particle-hole symmetry is broken. Thus, fluctuations in 
the amplitude of the gap cause charge density fluctua- 
tions, which couple directly to the potential $, see also 
Section IV for a discussion of this point. 

Our discussion below will be within the imaginary time 
Matsubara formalism. Here we will not address the im- 
portant point of relaxation mechanisms like spin-flip, 
electron-electron, and electron-phonon scattering that 
could be accounted for. We will simply assume that some 
relaxation mechanism is available which brings our super- 
conducting system into equilibrium with a big reservoir. 
We would like to emphasize that this assumption is by no 
means in contradiction with the main goal of our paper: 
to provide a convenient approach for studying dynami- 
cal and nonequilibrium phenomena in superconductors. 
Rather, it restricts the scope of the phenomena which 
can be effectively described with our methods. 

In what physical situations is our imaginary time for- 
malism meant to work? One such situation is quite stan- 
dard: a superconductor only slightly driven out of equi- 
librium, so that one can describe nonequilibrium effects 



within a linear response theory and express the results in 
terms of equilibrium correlation functions, which - after 
a proper analytic continuation - will describe the dynam- 
ics of the system in real time. Another important class of 
phenomena is related to (quantum) fluctuations of the su- 
perconducting order parameter, both of its modulus and 
its phase. Such fluctuations may - and in general do - 
involve virtual states with the electron subsystem driven 
far from equilibrium. Examples, are quantum phase slips 
in thin superconducting wiresEa and quantum tunneling 
of vorticeali. All such processes are also conveniently 
described within the formalism developed below. 

The analysis of the real time dynamics of a supercon- 
ductor with strong deviations of the quasiparticle dis- 
tribution function from equilibrium in general requires 
methods based on the Keldysh technique that teep track 
explicitly of the distribution function, see Ref.EIJ. Strong 
non-equilibrium real time dynamics is beyond the scope 
of the present paper and we postpone the corresponding 
discussion to a forthcoming publications. Here we would 
only like to point out that many features of our imagi- 
nary time analysis can be directly generalized and used 
also within the real time Keldysh technique. Thus, also 
for the real time case a lot can be learned already from 
the present imaginary time formulation of the problem. 

The calculation that leads to the effective actions 
Eqs. will be presented in the next Section. Its physi- 
cal content is discussed in Sections III (normal state) and 
IV (superconducting state). Where possible, we give the 
explicit forms of the propagators of the various fields in 
the hydrodynamic limit (in Sections III and IV), whereas 
the more general (and more complicated) expressions are 
deferred to Appendix B, where also the calculation of po- 
larization bubbles is outlined. The derivation of the Ward 
identities is given in Appendix A. n ,— ■ 

Parts of the present paper were implicit .ip-RefsxrE^I 
and to a somewhat lesser extent in Refs.OEJ IEj&O. 



II. MODEL AND DERIVATION OF EFFECTIVE 
ACTION 

The starting point for our analysis is a model Hamil- 
tonian that includes a short range attractive weak cou- 
pling BCS and a long range repulsive Coulomb interac- 
tion. We represent the latter in terms of the fluctuating 
gauge fields of electro-magnetism, V and A. The idea 
is to integrate out the electronic degrees of freedom on 
the level of the partition function, leaving-u s with .an ef- 
fective theory in terms of collective fieldsEiH^EMJ. The 
partition function Z is conveniently expressed as a path 
integral over the anticommuting electronic fields ip, ip 
and the commuting gauge fields V and A, together with 
a gauge condition. The Euclidean action reads 



dx 
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Here £(V) = — V 2 /2m — // describes a single conduction 
band with quadratic dispersion, A is the BCS coupling 
constant, a =f, \ is the spin index, and erii denotes the 
background charge density of the ions. In our notation 
dx denotes d 3 ~x.dr and we use units in which h and ks 
are set equal to unity. The field strengths are functions 
of the gauge fields through E = — W + (l/c)d T A and 
B = V x A in the usual way for the imaginary time 
formulation. 

We use a Hubbard-Stratonovich transformation to de- 
couple the BCS interaction term and to introduce the 
superconducting energy gap A = | Aje** 3 as an order pa- 
rameter 



exp ( A / dxip^ipiipiip^ ) = 
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V 2^ e -Jdx(\ 1 |A| 2 +Ai/> T '/u+Ai/>i</>T) (4) 



where the first factor is for normalization and will not 
be important in the following. As a result, the partition 
function now reads 

Z= [ X> 2 AWX> 3 AX> 2 *exp (-S - [ dx^Q^ 1 ^ , 



S [V,A,A} = J dx 
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where the prime on the integral denotes the restriction to 
a certain gauge choice for the electromagnetic potentials 
V and A. Below in Section III, we will sometimes use the 
Coulomb gauge V • A = in which the vector potential is 
completely transverse. In Eq. (JsJ) we have also introduced 
the Nambu spinor notation for the electronic fields 



and the matrix Green's function in Nambu space 

G - 1= ( d T -ieV + t(V~fA) A 
\A d T + ieV - £(V + fA) 

6 = 1 ° - 
1 F G 



(6) 



with normal and anomalous Green's functions denoted 
by G and F. 

After a final Gaussian integration over the electronic 
degrees of freedom, we are left with the effective action 



S eff = -Trln^- 1 + S* [^,A,A] 



(8) 



Here, the trace "Tr" denotes both a matrix trace in 
Nambu space and a trace over internal coordinates or 
momenta and frequencies. In the following "tr" is used 
to denote a trace over internal coordinates only. 



A. The Equations of Motion 

The Euler-Lagrange equations obtained by varying the 
action Eq. (||) with respect to V and A yield the two 
Maxwell equations that describe Thomas- Fermi and Lon- 
don screening, respectively. They read 



V • E = 4irie[n e — m] , 
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(9) 



Note that the ionic background contributes only to the 
charge density, and not to the current, if we describe the 
system in the frame where the ions are at rest. Both the 
electronic density n e and current density j e are expressed 
through the diagonal elements G and G of the matrix (in 
Nambu space) electron Green's functions Q. Explicitly, 



n e (x) = Tr[0cr 3 ] = G(x,x) - G(x,x) , 

j e (a;) = -Tr[(iVl + -Aa 3 )G] 
m c 



e 
rn 



(iV + -A)G(x, y) + (iV - -A)G(x, y) 

c c 



(10) 



y=x 



The matrices 01,2,3 are the Pauli matrices and below we 
will also use a± = \{p\ ± io^). 

The electronic density n e is a function of the chemical 
potential [i, and in the presence of particle-hole asym- 
metry also of the energy gap A. At zero temperature 
it satisfies n e (/x = e_F,A = 0) = rii. In general the 
electronic density can be expanded as n e (fj, + ieV, A) = 

n t + 2ieN Q V + 2N TA 2 H , with the density of states 

per spin Nq and the particle-hol&,asymmetry parameter 
T = d £ N(e F )/(2AA 2 ), see RefM The requirement of 
overall charge neutrality makes the electrostatic potential 
V a function of the energy gap, Va = — iTA 2 /e. Lon- 
gitudinal electric fields and deviations of the electronic 
density n e from the ionic density rij are screened on the 
Thomas- Fermi length scale \^ 2 F = 87re 2 7Vo. In the su- 
perconducting state, in addition the magnetic field B is 
screened on the scale of the London penetration depth 
A^ 2 = 47re 2 n s /mc 2 , where n s denotes the superfluid den- 
sity. 

Varying the action Eq. (||) with respect to A yields the 
BCS gap-equation for A 



A(at) = Tr[gcr_] = XF(x,x) 



(11) 



with the anomalous Green's function F. The gap equa- 
tion has a constant solution Ao = | Ao| exp(iipo) as well as 
more complex time and space dependent solutions, such 
as vortices. 



B. Perturbation Expansion 

The effective action in Eq. (^) is the starting point for 
an expansion around the constant saddle point solution 
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A = Aq, V = Va, and A = 0. We absorb the constant 
Va in the chemical potential fi from now on, so that it 
doesn't appear explicitly in the following. 

There are two ways of organizing the perturbation ex- 
pansion. In this section we will expand in V, A, and 
Ai = A — Aq, and to this end split the inverse Green's 
function in Eq. (J71) into an unperturbed part Qq 1 and a 
perturbation , according to 

00 = ( ? a ) ' 0fl = ( \ L -K- l ) ' 

jr^(l)V-ie7, l = zLL {V ,A}, (12) 
2 Vmc/ 2 toc 

where {., .} denotes an anti-commutator. In the following 
it is understood that the unperturbed Green's function 
has a chemical potential (i + ieV a and an energy gap A . 
Without loss of generality, we choose Ao to be real. 

At first sight the splitting we have just made seems 
not convenient and restricts severely the generality of our 
analysis, since we expand around a state with constant 
phase of the order parameter. However, this does not 
mean that we exclude, e.g., current carrying states for 
which the phase of A depends on coordinates in an es- 
sential way and is not small everywhere in the supercon- 
ductor. We show below that the actual parameters of our 
expansion are the gauge invariant linear combinations of 
the electromagnetic potentials and the phase of the order 
parameter, $ = V— ip/2e &ndv s = 2^(V</?— —A). Only 
these parameters (and not V, A, and ip separately) are 
required to be small within the framework of our analy- 
sis. Thus, also states that carry a current which is not 
necessary small can be described. This is a direct con- 
sequence of gauge invariance which plays an important 
role in our consideration. 

The other way of organizing the expansion is com- 
mented upon in subsection E. It involves a unitary gauge 
transformation of the fields, after which one expands di- 
rectly in the gauge invariant fields $ and v s , and a man- 
ifestly real perturbation Ai = |Ai|. With the help of a 
Ward-identity, the two expansions can be shown to be 
fully equivalent. For pedagogical purposes we postpone 
the corresponding discussion (see subsection E) and pro- 
ceed with the expansion. 

The trace of the inverse Green's function can be ex- 
panded in Q^ 1 using 

Trln^- 1 = Trln^ 1 + Tr ]T ^— ^ (SoSf 1 )" , (13) 

n=l 

and only terms of order n = 1 and 2 will be needed here. 

The n = 1 term in the effective action is Si = 
— ^'(GoGi 1 )- The explicit evaluation of the trace yields 

Si = -tr[K(G — G) - L(G + G) + A X F + AiF] (14) 



where we have used G(x, x) = — G(x, x) — —n e /2, and 
F(x,x) — Aq/X according to the gap equation. 

The second order contribution 5 2 = ^T^iGoQi 1 ) 2 
reads 

#2 = ^tr[GKGK + GKGK ~ 2FKFK + 2FLFL 

+GLGL + GLGL-2GKGL + 2GKGL-2FKFL 
+2FKFL + 2FAiGK - 2FA X GK + 2GAiFK 
-2GA1FK - 2FAiGL - 2FAiGL - 2GAiFL 

-2GA1FL + FA1FA1 + FA1FA1+ 2GA1GA1]. (15) 

As we are interested in contributions up to second or- 
der in the fields, it suffices to take K = —ieV in this 
expression. 

C. Longitudinal and Transverse Physics 

The next step is the evaluation of the traces in the 
expansion Eq. (|15|). For two Green's functions G and G', 
and two fields A and A' , the following identities hold: 

tr[GAG'A'} = J dqA(q)A'(-q){l} GGI ; 
tr[GAG'{W a ,A' a }} = 

= 2iJ dqA(q)A' a (-q){{p + q/2) a } GG , ; (16) 
tr[G{V a ,A a }G'{V b ,A' b }] = 
= -aJ dqA a (q)A' b (-q){(p + q/2) a (p + q/2) b } GG , , 

where a = x,y,z (repeated indices are summed over) 
and we have introduced the short hand notation q = 
(q, W/j), J dq = T J d 3 q/ (27r) 3 , as well as the bracket 
notation for polarization bubbles 

{B} GG , (q) = J dpBG(p + q)G'(p) . (17) 

Furthermore, we will split all fields in longitudinal and 
transverse components with respect to the momentum q 
by making use of the corresponding projection operators 

Pf = q a q b /c{ 2 , = 5 ab - q a q b /q 2 (18) 

that satisfy Pf = P u P L P T = P T Pl = 0, and P L +P T = 
1. With the help of the projection operators any vector 
field V a (q) can be decomposed into a longitudinal part 
V£(q) = P£ fc V b (q) = (V • q/q 2 )q a and a transverse part 
V£(q) = P% b V b \(i) = V a - (V • q/q 2 )<7 Q . Similarly, ten- 
sors T ab (q) are decomposed into T£ b (q) = Tr[P L T]P£ b 
and T£ b (q) = ±Tr[P T T]P£ b , where the extra half in the 
latter arises since TrPy = d — 1 = 2, in contrast to 
TrP L = 1. 

We use these considerations to simplify the polariza- 
tion bubbles that we encounter in the perturbation ex- 
pansion. By splitting (p + q/2) a into longitudinal and 
transverse components, we decompose the bubbles as 
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{(P + 9/2)a} 



GG' 



{Q/q 2 } 



GG 



,q 



{(p + q/2) a (p + q/2) b } GG , = {Q 2 /q 2 } GG , Pf + 



{(pxq) 2 /2q 2 } GG ,P^ 



(19) 



where in order to simplify notation, we have introduced 
Q = Q ' (P + <j/2)< In the following, we need only specific 
polarization bubbles for which we introduce the notation 

<70 = {1} GG (<?) , 5i = {Q/q 2 } GG (9) , (20) 
.92 = {Q 2 /q 4 } GG (?) , 93 = {(P x q) 2 /2q 4 } GG (9) . 

Analogously the fi,hi,ki denote {..}ff, {■■}gg> an d 
{■■}fg respectively. j-. 

Finally, as pointed out in Ref.E-3, it is advantageous to 
also split the fluctuations of the energy gap in real longi- 
tudinal and transverse components Ai = At +iA?. We 
will see in the next subsection that At and At appear 
in the effective action in rather different ways, related to 
their different physical nature. 



of combinations of the gauge fields and the transverse 
gap fluctuations as v s = ^[(VA T )/A — —A] and 

$ = y-^A T /A . 

The terms in Eq. ( pl| ) that couple the transverse phase- 
like gap fluctuations to the other fields assume a diagonal 
form with the use of the Ward-identity Eq. (|A9|). Little 
algebra shows explicitly that 

[A" 1 + h - f ]A 2 T + Ak Q V(q)A T (-q) 
-(4ze/m C )fc iq • A(q)A T (-q) = (4 e 2 A fc /^)[$ 2 - V 2 } 



+4mAohPf [{e 2 A a A h l 



b I 2 2\ a 61 

'm c — v.v.l 



which we use to eliminate At from the action. Further- 
more, the field strengths are invariant under gauge trans- 
formations, so that they may be expressed in $ and v s 
as 



|E| 2 =q 2 |<%)| ; 



2 2 



v«(?)| S 



[$(g)q • v s (-q) + *(-?)q • v a (g)] 



IBI 



D. The Effective Action 

The effective action up to second order in the fields At, 
At, V, and A is found by gathering the terms from So, 
Si, and S 2 (Eqs. (|5|, [l^ , |l5| ) ) . We split S e ff into a constant 
mean field part and a Gaussian fluctuation part S 2 ff . 
No first order contribution is present, as the terms ieriiV 
and —ien e V cancel by charge neutrality. We find 



S e ff = -Trlna o - 1 [Ao] + /3VA 2 /A, 
rE 2 + B 2 



(21) 



D off 



tfq 



A 2 T 



+ ( A L A q 



eV ^1 

mc 



A 



nm / e 




A 2 


2 \mcJ 




A L X 






A T 






eV 


: 


) 


eA b /mc J 


-Q 





where V denotes the volume of the system, (3 the inverse 
temperature, and we have introduced a matrix notation 
with 



M q = 



( h + f -iho -2iko 

iho h - fo -2k 

-2ik 2k fo - go 

\ -2q a h -2iq a ki iq a g 1 ? 

m ab = q 2 [(.93 + h)P T b + (.92 + h)Pf] 



2q*k l \ 
) 



2iq b k 



(22) 



In the above expression for the matrix M. q it is under- 
stood that all kernels are taken at momentum and fre- 
quency q. 

The physical content of the effective action Eq. ( |2l| ) 
can be brought out more clearly by "diagonalizing the 
matrix", i.e., rewriting Eq. ( ^lf ) in terms of the eigen- 
modes. To this end we introduce the superfluid veloc- 
ity and the chemical potential for Cooper pairs in terms 



-q 2 Pfv^(q)v b s (-q) , 



(23) 



which allows us to rewrite the terms related to A and 

V in terms of E a nd B. Finally, due to the Ward- 
identities Eqs. ( A7, A8|) the remaining terms in A and 

V are seen to vanish. Together with the propagator 
XA (q) = [2X- 1 + h (q) + h Q (-q) + 2f (q)} for the lon- 
gitudinal gap fluctuations and introducing the London 
and Joscphson susceptibilities Xt(?) = —8mAoki(q) and 
Xj{q) = 8e 2 Aofco(g)/icL' /J , we obtain the "normal" and 
the "superconducting" contributions to the effective ac- 
tion 

S sc [A L ,$,v s ] = -Trha^-^Ao^otO]) +/3VA 2 /A + 

+\ J <?9 (xaA 2 l + x,/$ 2 + Xlv 2 ) , (24) 
-Trln^- 1 [0] + 

XmB 2 ) . 



fE,Bl 



-- / dq { X eV z 



(25) 



Here we have also introduced the electric and magnetic 
susceptibilities Xe(q) — 2e 2 gi(g , )/(rrou; At ) and Xm(9) = 
2e 2 [ ff2 ( 9 ) + / 2 (9) - 9M - /3(9)]/(m 2 c 2 ). 

The terms in Eq. ( |2l|) that couple the longitudinal gap 
fluctuations to the other fields are nonvanishing only if 
particle-hole symme try is broken. Again by virtue of the 
Ward- identity Eq. (A9), they combine into the action 
S p h that describes the effects of particle-hole symmetry 
breaking 



Sph — 



aq(- X r(q)A L (q)A L (-q) 



(26) 



+ xkqWq) A U-q) + Xr (<?)q • v s (q)A L {-q) 



where we have introduced the susceptibilities xr(q) — 

ho(q) - h (-q), x|(g) = -2ie[k (q) + k (-q)], and 
X ^q) = 2[k 1 (q)-k 1 (-q)]. 
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Defining 

S* om [E,B] = J dq 
we obtain the final result 



E 2 +B 2 

8^ 



(27) 



S eS = S SC [A L , $, v s ] + 5 nm [E, B] 
+ S* ph [A L ,$,v s ] + S om [E,B] • 



(28) 

Using the standard definitions e = 1 + and /i" 1 = 

1 — 4ttxm, we arrive at the form Eq. (Mj quoted in the 
introduction. The expressions Eq. (p3-Eq) represent the 
main result of this Section, and a convenient starting 
point for any study of dynamical processes. 

Note that in arriving at Eq. fl28|) the terms ■izien i / e V 
from So and S\ have cancelled, since on the average the 



cancel. This point 
" 3 ^, where no coupling 



electronic and ionic charge 
was not appreciated in Rcfs. 
to electromagnetism was included and only the term from 
Si was found. 

We have decomposed the action into four parts: the 
superconducting contribution S sc , the normal metallic 
contribution S nm , the particle-hole symmetry breaking 
action S p h, and the action of the free electromagnetic 
fields S eia . Let us emphasize that the possibility of such a 
decomposition is a direct consequence of the Ward iden- 
tities. As these identities follow from gauge invariance 
only and do not depend on the presence and concentra- 
tion of impurities in a superconductor, we conclude that 
the splitting of the full action into four parts in Eq. ( psj ) 
holds not only for clean superconductors but rather for 
an arbitrary concentration of impurities. 



E. Gauge Invariance 

At this stage it is appropriate to discuss the con- 
sequences of gauge invariance for the effective action 
Eq. (IH) a bit more deeply, see also Appendix A. Inspect- 
ing Eq. (p8|), we observe that the transverse component of 
the energy gap At has completely disappeared from the 
effective action. This is just the Anderson-Higgs mecha- 
nismt3: the Goldstone mode At is "gauged away" and 
appears only within the combinations $ and v s . Since 
the electromagnetic field strengths E and B have been ex- 
pressed in terms of $ and v s , the integral over the field 
At in the partition function factorizes and contributes 
an irrelevant constant. 

The partition function Z can be represented in two 
equivalent ways. In the first, the Goldstone mode is ex- 
plicitly present and the four gauge field components are 
restricted by a gauge condition. In the second, the Gold- 
stone mode is "eaten" by the gauge condition and the 
four gauge field components are unrestricted. Explicitly 

J VA L VA T VVV 3 AeM~S cS {A L ,A T ,V,A}) = 

= y"pA L P$2? 3 v s exp(-S* cff [A L ,$,v s ]), (29) 



where the prime on the first integral denotes that it is 
supplemented by a gauge condition. In both cases 5 dy- 
namical degrees of freedom are present. 

We now return to the point raised in subsection B 
concerning the two possible ways of organizing the ex- 
pansion. Observe that the trace of the inverse Green's 
function, which is the starting point of the perturbation 
expansion, is invariant under unitary transformations 



Trln^T 1 = TilnUg^W 



(30) 



Choosing U(9) = exp(— ias9/2) with an arbitrary space 
and time dependent function 9{x), we obtain from the 
old Green's function Eq. (|?]), the new inverse Green's 
function Q~ x — exp(—ia30/2)Q~ 1 exp(ias9/2) that reads 



, i / d T — ie$ + £(V + imv s 



v e" ,( A\ 

= V e l9 A d T + ie$ - £(V - imv s ) ) ' ^ 

If A = |A|e* v , such a gauge transformation can be 
used to make the energy gap in the Green's function 
real by choosing 9 — tp. Instead, the superconducting 
phase appears in the chemical potential for Cooper pairs 
$ = V — tp/2e that replaces V and in the superfluid ve- 
locity v s = 2^r(V</5— ^fA) that replaces — eA/mc. Thus, 
we may identify <p with At/Aq from the previous sub- 
section. 

Note that the unitary operator U{9) is related to a 
local U(l) gauge transformation 



A -> e~ l6 A , (32) 



that leaves the field strengths invariant. Since the phase 
of all charged fields is rotated by a gauge transforma- 
tion, one also should replace the term ierii V in Eq. (§) 
by ierii&. The point being, that the ionic background 
charge density is eventually made up by particles and 
corresponding fields as well. 

After the gauge transformation with 9 = ip, the per- 
turbation expansion can also be done on the level of Q~ x ■, 
in terms of the real field Ai and 



(33) 



The whole derivation of the effective action is completely 
analogous and even slightly easier. In this way, one again 
recovers the result Eq. (|28|), now without making use of 
the assumption about small electromagnetic potentials 
and phase. These considerations emphasize the remark- 
able role of gauge invariance and conclude our derivation 
of the effective action for a BCS superconductor. 
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III. THE NORMAL METAL AND 
GINZBURG-LANDAU THEORY 

A. The normal metal 



As a first application of Eq. (|28|) , we will consider the 
normal metal limit for temperatures above the critical 
temperature Tq. If one puts Ao = in the normal metal 
and discards fluctuations of the energy gap, the electronic 
polarization terms can be expressed in terms of the lo- 
cally gauge invariant field strengths E and B only as is 
evident from Eq. (^) since in the normal metal S sc = 0. 
This contrasts to the superconducting case where also 
terms in <& and v s survive. For the normal metal we 
obtain 



5em — 



E 2 +B 2 



8tt 



rffoE 2 



(34) 



This action describes standard metal physics in the REA 
approximation, as is discussed for instance in Ref.EZl. 
Its analysis is most conveniently done in the Coulomb 
gauge. After an analytic continuation to real frequencies 
iuj^ — ► lu + i5 and |cj m | — > —iuj, the zeroes of the prop- 
agators describe modes of the electronic system. Two 
useful limits are the clean and dirty limit, in which 
0^,-Dq 2 » r" 1 and to^.Dq 2 <C t^ 1 respectively. Here 
D = vpr r /3 is the diffusion constant with a single par- 
ticle relaxation time due to impurity scattering r r and 
Fermi velocity vp. For frequencies of the order of the 
Fermi energy or the plasma frequency metallic systems 
are always in the clean limit. 

In the dirty limit, go is given by Eq. ( B12| ) and the part 
of the action related to the longitudinal electric field has 
the form 



Slim "I" S c 



(35) 

- l \u M \ + Dq 2 J 8tt 



kj,p, with Thomas-Fermi wavevector 



Since 87re 2 7V 

kxF, the low frequency part <C Dq 2 of this action 
describes metallic screening. In the opposite high fre- 
quency limit the second term of the action describes dis- 
sipation!^]: in terms of the conductivity a = 2e 2 NoD, the 
second part reads a~E 2 /2\u^\. 

In the clean high frequency limit ui^ 3> t~ 1 ,vf<1, 
the kernel go is given by Eq. (B7), and using go 
— N v F q 2 /(3uj 2 ), the longitudinal plasmon at frequency 
lu 2 = 4ire 2 n/m is recovered. 

The remaining part of Eq. ([mJ), related to the mag- 
netic and transverse electric fields, describes purely trans- 
verse physics. The kernel 33 can be approximated for 
low momenta and frequencies as 173 = (pp/3q 2 )go + 
^0/12 + • • •, and in the normal state the kernel 02 can 
be expressed through the Ward- identities Eq. ( A7 , A8 ) as 



92 = — nm /2q 2 + (moj / q 2 ) 2 go- In the dirty limit the bub- 
bles <72 and <?3 almost cancel. The remaining action for 
the transverse vector potential in the dirty limit is 



S'n 



(36) 



\uj^\ + Dq 2 ) 8tt 



The low frequency limit describes the normal skin-effect, 
i.e., iu> — c 2 q 2 /(47rcr), for wavelengths larger than the 
mean free path I = vft t . 

In the opposite clean limit and for small frequen- 
cies lu^ <C vpQ, we find explicitly 53 — 32 ~ No/ 12 + 



TTNolto^/iAvFq 3 ) 



The dispersion is now different, 



iio = 4upA 2 (0)|(?| /37T, with the zero temperature Lon- 
don length given by A^ 2 (0) = Aire 2 n/mc 2 , and is related 
to the anomalous skin-effect. Landau dia-magnetism is 
recovered from the small constant term in (73—32 = 
Nq/12+- ■ •, and this is a (vp/c) 2 correction to the "l/87r" 
in S em in Eq. (§7|): fx' 1 = 1 + l/(4A|(0)p|). Pauli para- 
magnetism is not present in the action Eq. (^) , since we 
did not include a Zeeman coupling in the original model 
Eq. (^). Finally, two transverse light modes with disper- 
sion oj 2 = ui 2 + c 2 q 2 are present in the high frequency 
clean limit. 



B. The Ginzburg-Landau Expansion 

In the normal metal close to Tc , fluctuations of the or- 
der parameter A can be studied using Ginzburg-Landau 
theati'^lcl or its time dependent (TDGL) generaliza- 
tionBEl "El. Within our present formalism, the TDGL 
effective action is readily derived starting from the ac- 
tion Eq. ((H). The expansion in $ and v s is already 
performed and we need only to expand all terms in pow- 
ers of A. In addition, we need parts of the third and 
fourth order terms in the expansion of the inverse Green's 
function Eq. (p"3|). They are calculated using the nor- 
mal metal Green's functions with zero energy gap, see 
App. B3. Th is procedure is quite standarde2l (see also 
Eqs. (B7, B12| )). The superconducting part of the action 
takes the form 



8T 



A| 



(37) 



+ N Q J da (A [in -£ 2 (0)V 2 +4m 2 £ 2 (0)v 2 

b, 



-T(d T - i2eV) + 2e 2 6$ 2 



A+2IAI 4 



Here we have introduced the coherence length £ 2 (0) = 
-fxD in the dirty limit and £ 2 (0) = (7C(3)/487r 2 )(^/T) 2 



8T 



in the clean limit. The coefficient b = 7C(3)/(8tt 2 T 2 ) is 
derived in Eq. ( |f31~6|) . The small coefficient T = N^/2\N 2 



arises due to particle hole asymmetrycil and is usually 
neglected. The term containing this coefficient describes 
the small difference between electronic and ionic densities 
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resulting from the fluctuations of the order parameter. It 
occurs in the quoted gauge invariant form, since the time 
derivative from the Cooperon combines with the term 
~ A 2 V from the third order expansion. We will see in 
the next Section that the same term arises also at low 
temperatures. 

The term involving the second order space derivative 
can be combined with the v 2 A 2 term from the fourth 
order terms in the expansion into one gauge invariant 
second order derivative £ 2 (0)|(V - i^A)A\ 2 . The $ 2 A 2 
term, however, does not straightforwardly combine with 
time derivatives into gauge invariant time derivatives. 
The reason is that close to Tc dissipative and Hamilto- 
nian frequency dependences mix. As an example, the dis- 
sipative |cj^| termla in Eq. (|37|), which turns into the dis- 
sipative time derivative of the real time TDGL equation 
after an analytic continuation, clearly cannot be made 
gauge invariant. Since the second and higher order time 
derivatives are usually irrelevant as compared to the dis- 
sipative jw^l term, we do not include them in Eq. (|37|). 

Let us also note that the expression Eq. |37]) is cor- 
rect only in the limit of low frequencies and wave vectors 
uif!, Dq 2 < 4irT. For larger frequencies the expression be- 
comes more complicated, as we can no longer expand the 
kernel ho in u^/lnT, e.g., #(1/2 + w m /4ttT) ~ *(V 2 ) -> 
■klo^/KT (fy is the digamma function, see App. B). Also 
the gradient terms in Eq. ( |37| ) should be modified in this 
case. As the corresponding expressions turn out to be 
quite tedious we do not present them here. For some 
problems, however, these modifications become signifi- 
cant, especially because the validity of the GL expansion 
Eq. ( |37| ) is restricted to temperatures T ^ Tc, in which 
case w M is never really smaller than 4-7rT. 

The action S nm is also important for the description 
of the dynamical properties of the superconductor and 
describes dissipation in the "sea" of the remaining normal 
electrons. It explicitly depends on the order parameter, 
since also the polarization bubbles can be expanded in 
A, giving rise to additional contribu tions. Expanding 
the bubble go + fo ( see Eqs. (B11.B11)) up to the second 
order in A, we Q b jt. ai ji in the limit of small frequencies 
and wave vectoral^'E^I 



AT Dq 2 Ar ^A 2 
9o +fo = -Nqt-^t - N - 



4K|T 



(38) 



For the part of the effective action concerned with the 
electric field, we find 

f * / 47R7 2ir 2 e 2 N A „\ E 2 

Snm + S cm = dq(l + — + „. ° A 2 — . (39) 



8tt 



The dependence of the action 5 nm on A is important 
in dirty superconductors and accounts for the Maki- 
Thomps on fl uctuation enhancement of the conductivity 
near TccJiiS. The Azlamazov-Larkin fluctuation correc- 
tion to the conductivity is already present in the TDGL 
actum, S sc Eq. ( |37| ) due to the presence of the v 2 A 2 
tcrmE-3. 



Extensions of the elegant TDGL action into-the super- 
conducting phase have turned out to be hardEEl. Only in 
thej-aresence of a large amount of paramagnetic impuri- 
tiese 2 ] or very close to TcO is this possible. 



IV. DYNAMICS AT LOWER TEMPERATURES 

A. Electromagnetism 

We will now evaluate the contents of Eq. ( |28| ) in the 
superconducting state with an equilibrium gap Ao. We 
first focus on the parts involving E and $ that are re- 
lated to the electric field and the Josephson relation. For 
the prefactor of <i> 2 , we take the bubble 4Aofco/itJ Al = 
Non s /n at zero frequency and momentum. At small fre- 
quency and momentum the bubble that multiplies E 2 is 

reads 



Nan n /(nq ), see Appendix B. The action 



S e s = J 3q ( 



_L + e2jV ° n " 

8tt 



r 



E 



e 2 7V — $ 2 
n 



(40) 



It describes metallic screening of the electrostatic po- 
tential V with the full electronic density n, to which 
both terms in Eq. ( [40| ) contribute, and superconducting 
screening with superfluid density n Sl only through the 
second term, to enforce the Josephson relation V = <p/2e. 
At higher frequencies and momenta, weight shifts from 
the term ~ n s to the term ~ n n in such a way that the 
plasma frequency and Thomas- Fermi screening length re- 
main constant. The higher order frequency and momen- 
tum dependence of the kernels is not easy to extract. Fol- 
low frequency and momentum, one typically finds correc- 
tions of order ui^/ A and vpq/A or Dq 2 /A . In the op- 
posite limit, for high frequencies and momenta, the ker- 
nels reduce to their normal state form, and we recover dis- 
sipation in the dirty limit for frequencies A( 



< 



< 



(cf. Eq. (p5|) ) . In the clean limit such an intermediate 
frequency regime does not exist. 

We now turn to the parts of Eq. (^8[) that are re- 
lated to the magnetic field and the superfluid velocity. 
To lowest order in the external momentum g 3 g 2 ~ 
(p 2 F /3q 2 )g + ■■■ and f a » f 2 » (p 2 F /3q 2 )f + so 
that the combination g% + fa — g 2 — fi is equal to zero 
in the q — * limit. What remains is the evaluation of 
— 4mAofci = mn s /2, and we obtain 



Sefi = 



— —v. 



(41) 



This action describes transverse screening of the mag- 
netic field in a superconductor and is related to the Lon- 
don theory. 

Summarizing: at high frequencies lo^ 3> Ao and mo- 
menta q 3> £ -1 the electromagnetic properties are those 
of a normal metal. At low frequencies and momenta a 
superconductor screens, in addition to electric fields, also 
magnetic fields. 



8 



B. Dynamics of the Energy Gap 

It is clear from Eq. that the dynamics of the fluc- 
tuations of the amplitude of the energy gap A l are gov- 
erned by the combination \A = [2A _1 + 2f (q) + h (q) + 
ho(— q)]. In the clean limit we obtain to lowest order in 
co^/Ag and vpq/Ao at zero temperature (see Appendix 
B) 



S sc = N Q I dq(l+±-^ 



2 2 

36A 2 



A! 



(42) 



This result is often quoted in the literature as evidence 
for the existence of a TDGL-like theory with a second 
order time derivative at zero temperatureEH. Its use is 
restricted, however. The mode that it describes has dis- 
persion lu 2 = 12AQ+(l/3)u|,q 2 , i.e., it is gapped to an en- 
ergy where the expansion in w M /Ao does not make sense. 
The correct mode of the amplitude of the energy gap is 
heavily overdamped due to the coupling to particle-hole 
pairSpr-aad starts off at frequencies 2A as discussed in 
RefsHEi For driven situations at small frequency and 
momenta on the other hand, the expansion in Eq. (G^) 
is useful. At nonzero temperatures, the exact frequency 
and momentum dependence of xa is rather involved un- 
fortunately. 



C. Particle- hole asymmetry 

We now turn to the particle-holc symmetry break- 
ing action Eq. (|2^). To lowest order in the frequency 
and momentum, we obtain Xv(l) = ho(q) — ho(~q) = 
2NdTw lt , xh<l) = -2*e[Ao(g) + k (-q)] = -4ieN TA , 
and Xrio) = 2[fci(<?) - h(-q)] = 0. Together with the 
term |Ai| 2 V from the third order expansion of the effec- 
tive action, we find 



S ph = -2ieN T J dx<5> (| A + Ai| 2 - |A | 



(43) 



as announced in the introduction. Apparently, this con- 
tribution to the action is independent of temperature and 
the mean free path. Close to Tq exactly the same term, 
with Ao = appeared already in the TDGL expansion 
(see Eq. (|37|)). The coupling constant V of this term is 
usually small, and S p h is irrelevant, except for inhomo- 
geneous problems related to vortex motiontl. 

In the core of a vortex the energy gap goes to zero, 
and this local variation of the gap induces a local charge 
density modulation, for which the action S* p h contains the 
source termE3. Moreover, due to the singular phase field 
around a vortex, together with the suppression of the gap 
in the core, the action S p h gives rise to a small additional 
force per length on vortices, — 27riVorA 2 |V£ X z, which is 
proportional and perpendicular to the vortex velocity 
(z is the unit vector along a vortex in the direction of 



the magnetic fieldJiiS. We do not find evidpice for the 
much larger force ~Tin s 'V[ J x z found in Ref.efl. 

Note that this "topological" force is only one contribu- 
tion to the several different forces on a vortex; for a pa- 
view of the other forces on a vortex- line, see, e.g., Rcf.c3. 



D. The Uncharged Limit 

The limit where the electronic charge vanishes has re- 
ceived some attention recentlyEj~t3. Although this case 
is realized in superfluid 3 Hc, the different order parame- 
ter symmetry makes any s-wave considerations less use- 
ful. Furthermore, the interactions between uncharged 
3 He atoms is very different from the electron-electron in- 
teractions as described by our starting point Eq. (||). In 
particular, 3 He atoms are neutral and in the 3 He system 
no background charge is present. 

For completeness, however, we also discuss the un- 
charged limit of our model Eq. (|^). Putting e — > in 
Eqs. (|2jj2^), we find to lowest order in momentum and 
frequency for the phase part of the action 



S sc + Sph = N / dq ( - iTAlip 



4?? 



9" 



f (V^) 2 



(44) 



which gives the standard acoustic Bogoliubov-Anderson 
mode with velocity v — vp/\/3. The velocity is tem- 
perature independent, and non-critical even around Tq. 
Note that the phase actios, Eq. ( (f4| ) is different from the 
ones pbtained in Refs.E3~E3. In particular, in contrast to 
Refs.Ej Ea no large topological term in e (p/2 is present, 
only a much smaller term proportional to the particle- 
hole asymmetry. The difference can be traced to the near 
perfect cancellation between ionic and electronic charge 
densities in our case. 

For superconducting Bose-liquids a large topological 
term is present in the effective action, since the bosonic 
field itself can be taken as the order parameter at low 
temperatures and as a result the phase is dual to the den- 
sity. However, for Fermionic superconductors one rather 
expects that the phase of the order parameter is dual to 
the amplitude of the order parameter, i.e., to the energy 
gap. As a consequence, instead of a term in e ip, we expect 
a term proportional to iA 2 (fi in the effective action. This 
is just the content of Eqs. (^3|,^4|), which shows that the 
constant of proportionality is given by the particle-hole 
asymmetry parameter T. 



V. CONCLUSION 

We have reviewed the derivation of the effective the- 
ories for BCS superconductors and discussed the corre- 
sponding dynamics of electromagnetism and the ampli- 
tude of the energy gap. Our main result Eq. (^8j) is a good 
starting point for investigations of quantum dynamical 
and statistical problems in BCS superconductivity. 







We have stressed the role of gauge invariance and the 
corresponding Ward identities that express particle num- 
ber conservation. Although a perturbation expansion can 
violate gauge invariance, the Ward identities allowed us 
to obtain explicitly gauge invariant results. In particular, 
we have demonstrated how the Anderson-Higgs mecha- 
nism occurs within BCS theory. In contrast, .the role of 
Galilei invariance that was stressed in Refs.t3t3 does not 
seem to play an important role in real BCS materials. 

Furthermore, we included the effect of particle-hole 
asymmetry in our considerations. We find a small topo- 
logical term proportional to the particle-hole asymme- 
try, that leads to an additional Hall-force on vortices 
(apart from the Kopnin-Kravtsov, .Magnus, and Ior- 
danksii forcesc3) as discussed in RcfsJj'E2l. Also, we have 
seen that the structure of the theory is essentially the 
same in the clean and dirty limits. In particular, the 
prefactor of the topological term does not depend on the 
electronic mean free path. The main difference between 
the clean and dirty limits is the presence of an interme- 
diate dissipative regime A < o;„, vpq < t^ 1 in the dirty 
limit. This difference shows up for instance in the .quan- 
tum dynamics of low dimensional superconductors!! 
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APPENDIX A: THE WARD IDENTITY 

An important step is the derivation of the Ward iden- 
tity related to gauge invariance. On the level of vertex 
functions this Ward identity is discussed for instance in 
Rcf.El. Here we will derive the Ward identity on the level 
of the Green's function in order to obtain relations be- 
tween the polarization bubbles. Our derivation holds for 
an arbitrary concentration of impurities in a supercon- 
ductor. 

In the clean limit the Ward identity also holds as an 
algebraic identity. For the case of normal electrons, with 



single particle Green's function G{p) 



-£ p ] \ it is 



G[p) - G(p + q) = (i^ + Q/m)G(p)G(p + q) , (Al) 

where Q/to = £(p + q) — £(p). 

A general way of establishing the Ward identity is to 
consider the change in the Green's function upon rotating 
the electronic phase by ip. We have on the one hand (we 
expand in ip) 



G<j>(x, x) = e 



Is _ i<t>(x)a 3 /2 



G(x, x 1 ) 



-i<fr(x')a 3 /2 



so that 



SG(x,x') = ~[<f>(x)*3g(x,x') - g(x,x')* 3 <f>{x')} (A2) 



dqdpe 



iqx ip(x — x ) 



whereas on the other hand 

Q-\ x , x>) = e^^g-^x, x')e^ x '^/ 2 

= g- 1 (x,x') + 5g- 1 (x,x') 

so that 

SG~Hv, v') = 8(v -»')(- 5^3 + ^{v, v^}i + 

+i<f>Aa+ — i(j)A.a- 

and 



Sg(x,x') = - / dydy'g(x,y)Sg- 1 (y,y , )g(y',x') 



= — I dqdpe 



iqx ip(x — x ) 



Jqyp-\-q 



UL a3 - -rLl + i (Aa+-Aa.) 
2 2 to 



(A3) 

g v ■ 



Comparison of Eqs. (A2],A3j) leads to 



C3^p — g p + q &3 = 
g p+q [iui^a 3 + Q/mi — Aa + + Aer_] g p . (A4) 



The restriction to the normal metal yields again Eq. ( Al ) . 
In the superconducting case the Ward-identities are 

slightly more complicated; The upper left and upper 

right components of Eq. (|A4|) read explicitly 



G(p) - G( P + q) = (iw p + Q/m)G{p)G(p + q) + 
+(-icJ M + Q/m)F(p)F(p + q) - 
- 2AF(p)G(p + q) + 2AG(p)F(p + q) , 



F(p) + F(p + q) = (i 



Q/m)F(p)G( P + q) 



+(-iuj„ + Q/m)G(p)F{p + q) - 
-2AG(p)G(p + q) + 2AF(p)F{p + q) . 



(A5) 



(A6) 



These identities can be used to generate the Ward iden- 
tities for the electronic polarization bubbles, by tracing 
them over the internal momentum and freq uenc y p to- 
gether with some function. The trace of Eq. (A5) with 1 
and Q/m immediately gives 

= iuj^goiq) - f (q)} + q 2 /m[g l (q) + fi(q)] 

+2A [k (q) - k (-q)] (A7) 
-n/2 = iw„[ffi(?) - fi{q)} + q 2 /m[g 2 (q) + f 2 (q)} 

+2A [k 1 (q) + kti-q)] , (A8) 



= g(x,x') + sg(x,x / ) 



and the trace of Eq. ( A6) with 1 yields 
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A A 1 = A [/o(g) - h (q)] + iuj^k (-q) 
-(q 2 /m)ki(-q) . 



(A9) 



The last three Eqs. (|A7, Aq, A9|) are used in the text in 



Section 2. They are a result of gauge-invariance (parti- 
cle number conservation) and hold also after an impurity 
averaging procedure. The identities can be simplified fur- 
ther using fi = which holds by symmetry. 

The previous discussion was based on the gauge sym- 
metry that is generated by exp(i<73</>/2). Along com- 
pletely analogous lines, the invariance with respect to 
rotations by exp(il</>/2) leads to 



i^nfo(q) = Ao[fc (g) - k (-q)] > 
-(q 2 /m)f 2 (q) = A [k 1 (q)+k 1 (-«)]. 



(A10) 



These last two identities are not imp ortant in the deriva- 
tion of the effective action Eq. (f2q). They do, however, 
reduce the amount of work needed to explicitly evaluate 
all the different polarization bubbles. 



For the other bubbles we obtain similar expressions with 



s g = 


[(E' + 




EE') + 


iu^'E - £E')\ , 


N g = 


[(E>- 


+ 


EE') - 


iu^'E + ZE')] , 


Sk = 


l(E' + 


E)£A + 


iu^EA] 




N k = 


[(£'- 


E)tA- 


iu^EA] 




s h = 


HE' 


+ E)(tf' 


+ EE') 


+ E + ZE*)] , 


N h = 


HE' 




-EE') 


-iu^'E-ZE')]. (B3) 



The remaining integral over momenta can in general not 
be given in closed form. Let us therefore consider the 
simple limits. For external momentum and frequency 
much smaller than the energy gap and for u>^ <§; vpq, the 
kernel fo reduces to 



fo = N Q J dZ 

_ N n s 
2 n 



\ 2 



■2E 2 \2E 



1 



[1 - 2f E ] + d E f(E) + 



(B4) 



For the others we obtain 



APPENDIX B: THE POLARIZATION BUBBLES 

In this Appendix the polarization bubbles that are 
used in the text are discussed and results in several limits 
are summarized. 

1. Clean Limit 

We evaluate the kernels by doing the sum over Mat- 



subara frequencies first. The notation E = \/£ P + Aq 



and E r = y J C P + q + ^ s use d, as well as J cfcl to denote 
a normalized angular integration. In our notation the 
unperturbed Green's function in momentum space reads 
explicitly 



G F 
F G 



-iu v + £ p A 
£ 2 + A 2 V A o - - £ P 



(Bl) 



For the bubble /q we obtain by standard contour inte- 
gration and ordering terms 

/(i 3 p 1 v 
J2^)3 p F ( P + q ' Wv + W M) i? (P) Uv) 



1 ( [l-fE>- J E \ 

2EE'\uj 2 + (E' + E) 2 



Nf) . (B2) 



[fw - [e\ 
-tf + iE'-E)*- 



Here fp = f{E) is the Fermi function, and Sf and Nf 



arc 



S f = (E' + E)A 2 ; Nf = (E' - E)A 2 



9o 



2 n 



Nq-ti- n r , 



1 



A 

n ta 



~2 n 

Np n s iu^ 

4 n A 



ho = — — + — h NoTiuj^ 

k 



(B5) 



The first order expansion of the kernels in frequency and 
momentum at zero temperature reads 



JO 7, 



90 = r- 



No 
2 

No 
2 



Np 
12A 2 
Np 
12A2 



V F 2 

— q 

3 H 



h = -- + — + NoTiLUf, 



k = N TA + — -r 11 

4 A 



Nq 
6A 2 
1 

6Af 



(B6) 



In the normal metal T > Tq and A = the expressions 
simplify considerably. For q <§C 2k p we have 



fo = k = 
90 = -N (l 

w<c VFq ■ s 
w> VFq ■ p 



(B7) 




\VFqJ 



Np ( VFqV 

3 V w m / 



h = No 



7T |CJ M | + 



7C(3)^| „2 



T 



+ »Tw M - In 
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The bubbles are related to one another by the exact 
Ward identities, and also by approximate identities that 
are good in the limit of low external momenta. As an 
example we discuss the relation between 173, gi, and go- 



,93 = {(pxq) 2 /2<? 4 } GG 



52-{[q-( P + q/2)] 2 /9 4 } GG « 



^ q Vf {1} 



GG 



3q 2 



.9n 



(B8) 



Here we have used (1/2) J <i(9 sin cos 2 9 = 1/3 and 
(1/2) f£ d9 sin sin 2 9 = 2/3 as well as the fact that inter- 
nal momenta can be taken at the Fermi energy. In doing 
so, one makes a slight error and it can be shown that the 
leading term in the difference (73 — 32 is Nq/12, which 
is responsible for Landau diamagnetism as discussed in 
Section III. 



2. Dirty Limit 

We will use the notation oj = lu v , lu' = lu v + lu^, W — 
+ A 2 ,, and W — \Jlu' 2 + A 2 ,. In the presence of 
impurities we replace all frequencies and the gap by lu — 
t\lo and A = ?7Ao, with 77 = [1 + l/(2r r W)], in the 
singles-particle Green's functions. In particular W = W + 
1/2tM. 

Furthermore, we will restrict ourselves to the limit q <C 
2k F and we put Q = xv F q, with x — cos 6?. Although 
formally one should first sum over the internal frequency, 
in the dirty limit it is more convenient to integrate over 
energy first. By reversing the order one only misses a 
constant —No in the expression for go, which is added 
by hand later on. The integral over energy £ p by contour 
integration is straightforward and the angular integration 
gives an arctan. We find for fo 



ladders is obtained from self-consistencylll by s impl y 
dropping the r^T 1 from the right hand side of Eq. ( BIO ). 
We find the final expression for the disorder averaged 
polarization bubble 



f°=« N °pY, ww , (w 



A 2 
^0 



W + Dq 2 ) 



and analogously in the same limit Dq 2 ,o; A1 Ao 
obtain for the other kernels 



we 



.9o = -N 



WW - LUUj' 



=4e 



W + Dq 2 ) ' 
-N (WW + lulu') + NqWW'ilu^ 



WW'(W + W' + Dq 2 







E 



No(-iLuA ) 



NqWW'A 



WW'(W + W + Dq 2 ) 



(Bll) 



A considerable simplification occurs for T > Tc when 
Ao = 0. In this case the remaining sums over the internal 
frequencies can be carried out and yield differences of 
digamma functions. The low momentum and frequency 
expansion gives 



fo = k = 
5o = -N c 



(B12) 



Dq 2 



Wu n 



h = N 



7T \luJ +Dq 2 



T 



iTLU a — In 



2e 1 uu D 
ttT 



In this limit, the kernels go and ho are nothing but the 
Diffuson and Cooperon. For temperatures T < Tc no 
simple expressions are available. However, the bubbles 
at zero external momentum and frequency are known 



/o(Q) 



2 p ~ (w 2 



A 2 
^0 



A 2 )3/ 2 



No n s 
2 n 



(B13) 



See Eq. (B5) for the other kernels. 



fo = / dp- 



A 



An 



ul + ^ + A 2 (^+^) 2 +£ 2 +q +A 2 



(B9) 



A 2 
^0 



(3 ^ v F qWW 



arctan 



v F q 



W + W + T r r 1 



In the dirty limit when AoT r <C 1, the arctan function 
can now be expanded in v F qr r — ql and Lu^r r , and to 
leading order 



1 , ( v F q 

arctan =- 

v F q \W + W + T r r 1 

1 



W + W + Dq 2 + Tr 



(BIO) 



where D = T r v F /3 is the diffusion constant. The full 
bubble including the vertex correction due to impurity 



3. Higher Order Bubbles 

In Section 3 we need the diagrams with 3 and 4 Green's 
functions at zero external momenta for the construction 
of the Ginzburg-Landau functional. The third order con- 
tribution Iq is independent of the mean free path 



dp[G(p)G(p)G(p) - G(p)G(p)G( P )} 

-2£ 



= _ 2iV 'iv— r dxx2 

: -N T , (B14) 



-N^hx 



2e<Lu D 
ttT 
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and is proportional to the particle-hole asymmetry of the 
problem, whereas the combination 1$ is 

I+ = J dp[G(p)G(p)G(p) + G(p)G(p)G(p)} 



. (B15) 



The fourth order contribution is 

h = J d P G 2 (p)G 2 (p) = N ± £ j #i H lp )2 

- n - _L_ f°° At dx 

- 7Vo / jtrKi 3 i_ 00 ^(i+ a; 2 ) 2 

Here we have introduced the usual Ginzburg-Landau pa- 
rameter b. Finally, below Tq we need the combination 

If = jd P [G(p) - G(p)][3F(p)F(p) + G(p)G] 

« -2iV r . (B17) 
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